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Dynamics of flexible flbers in shear flow 
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Dynamics of flexible non-Brownian fibers in shear flow at low-Reynolds-number are analyzed 
numerically for a wide range of the ratios A of the fiber bending force to the viscous drag force. 
Initially, the fibers are aligned with the flow, and later they move in the plane perpendicular to 
the flow vorticity. A surprisingly rich spectrum of different modes is observed when the value of 
A is systematically changed, with sharp transitions between coiled and straightening out modes, 
period-doubling bifurcations from periodic to migrating solutions, irregular dynamics and chaos. 


In the literature, there has been for decades a lot 
of interest in analyzing the dynamics of flexible non- 
Brownian and Brownian fibers in shear flow at low- 
Reynolds-number: experimentally, theoretically and nu¬ 
merically [IH 3 . Recently, migration and rheology related 
to Brownian motion in shear and Poiseuille flows have 
been extensively investigated [Ml. The growing inter¬ 
est in this field has been motivated by its importance for 
applications for polymers, proteins, DNA or other bio¬ 
logical systems. 

On the other hand, the coil-stretch transition of flex¬ 
ible polymers discussed by de Gennes 0 was followed 
by numerous fundamental theoretical and experimental 
studies of buckling instabilities of non-Brownian fibers in 
different ambient fluid flows [Tsl - I^ . The results demon¬ 
strate a complexity of the dynamics of flexible fibers, and 
its sensitivity not only to the fiber bending rigidity, but 
also to many other different factors, such as e. g. cur¬ 
vature of the flow or presence of walls [ 1 ^, very large 
length of fibers leading to topological changes of their 
shapes [llj, [ 2 ^ or a spectrum of initial conditions leading 
to various modes of the three-dimensional motion Q. 

Therefore, in this work we focus on a simple system: 
a single non-Brownian deformable fiber in an ambient 
shear flow of a fluid with viscosity ry. The goal is to 
investigate the dependence of the dynamics on the fiber 
flexibility. We try to minimize or eliminate other sources 
of a complex behavior. Therefore, the flow is unbounded 
(there are no walls), the fiber is moving freely, and its 
aspect ratio is not large. Initially, the fiber is straight, 
aligned with the flow, and there is no constraint forces 
acting on any of its segments, what leads to the fiber 
motion entirely within the plane y = 0 perpendicular to 
vorticity of the ambient shear flow with velocity vq = 
'fz'k. 

A fiber is modeled as a chain of N identical solid spher¬ 
ical beads of diameter d [2^ . The length, time, force and 
velocity units are chosen respectively as 

d, t = 1 / 7 , /o = 7rr7d^7, vc = dlT. ( 1 ) 

The centers of the consecutive beads are linked by 
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springs with the equilibrium length l^d so small that the 
consecutive beads almost touch each other. The spring 
constant kfo/d is so large that the fiber’s length practi¬ 
cally does not change. 

At the equilibrium, the fiber is straight. At a deformed 
configuration, there appear a bending force exerted on 
each bead, proportional to the (dimensional) bending pa¬ 
rameter Af^d? [3,[H,[23,0. Our goal is to analyze how 
the dynamics of fibers depends on their dimensionless 
bending stiffness A. 

We assume that the Reynolds number is much smaller 
than unity, and the fluid velocity and pressure satisfy 
the stationary Stokes equations [^. Velocities of the 
beads at given positions are evaluated with the use of 
the Hydromultipole numerical fortran code, based 
on the multipole method of solving the Stokes equations, 
with the multipole expansion truncated at an order L 
[ 2 ^. The time-dependent positions of the beads are 
determined by the adaptive fourth order Runge-Kutta 
method. In addition, the time-step is reduced by a fac¬ 
tor of two, if potential overlaps are detected, what allows 
to avoid spurious overlaps (26l[ . 

We choose the following values of the parameters. 


N = 10, /o = l-0I, fc = 2000, 0.3<A<50, L = 2. (2) 


While moving along the flow, the fiber center-of-mass 
oscillates across the flow, with or without a systematic 
lateral migration, depending on the bending stiffness A. 
The fiber tumbles, in a similar way as a rigid elongated 
body [ 13 . While it turns, its shape evolves accordingly. 

To determine the lateral motion, we evaluate the 

time-dependent z-component of the fiber center-of-mass 
position. Examples are shown in Fig. [TJ 

To describe the tumbling, we first define the end-to-end 
vector as the difference of the positions of the centers of 
the last and first beads. Then, we determine the flipping 
instants t/ as the moments when the end-to-end vector is 
oriented along 2 . At a flipping instant, the center-of-mass 
lateral position is denoted as 2 / = Zm(tf). The tumbling 
time r associated with tf is defined as the difference be¬ 
tween tf and the next consecutive flipping instant. 

To analyze deformation of the fiber shape, we evaluate 
the time-dependent fiber curvature K{t) and fractional 
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a) A=0.55 


b) A=0.7 


C) A=0.9 


d) A=1.14 



FIG. 1: Evolution of the fiber center-of-mass position Zm, curvature k and fractional compression a depends on the 
fiber stiffness A. Red dots indicate flipping times, and black circles - times of the corresponding snapshots in Fig. [21 


compression a{t), defined as 
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( 3 ) 


where is the radius of the circle determined by the 
centers of three consecutive beads, and 6{t) is the fiber 
end-to-end distance (i.e. the length of the end-to-end vec¬ 
tor), initially equal to the equilibrium value {N — 1)^0- 
Examples of K{t) and a{t), shown in Fig. [H and the cor¬ 
responding evolution of hber shapes, illustrated at the 
snapshots from the simulations in Fig. |2l provide exam¬ 
ples of different modes of the dynamics, which now will 
be introduced and discussed in details. 

We first observe that fibers in Fig. |2ja,b) all the time 
stay coiled, but fibers in Fig. |2]^c,d) straighten out from 
time to time. For example, the fiber with A = 1.14 
straightens out along the flow, then forms a compact 
S-shape, straightens out again, and so on, as shown in 
Movie (d). In contrast, the fiber with ^ = 0.7 always stays 
coiled - it does not straighten out, as shown in Movie (b). 
Its end-to-end vector rotates, with a well-defined tum¬ 
bling time, but the trajectories of all the fiber beads stay 
close to each other, and the motion of the beads resem¬ 
bles tank treading observed for vesicles [^ . 

The differences between the coiled and straightening 
out modes are illustrated in Fig. [TJ For the fibers which 
straighten out, the curvature and fractional compression 
oscillate between a large value (corresponding to a com¬ 
pact shape) and almost zero (corresponding to almost 
linear configuration). For fibers which do not straighten 
out and stay coiled, K{t) and a{t) also oscillate, but all 
the time their values remain large. 

To find the stiffness threshold for the coiled - straight- 
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FIG. 2: Shape evolution of fibers in shear flow for the 
fiber stiffness and tumbling time (a) A = 0.55, r = 370, 
(b) A = 0.7, t = 352, (c) A = 0.9, t = 2I9, (d) A = 1A4, 
r = 206. Snapshots are taken at equal time intervals. 
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ening out transition, we perform simulations for a wide 
range of values of A, and analyze the minima K^m and 
o^min of the hbei' curvature K(t) and fractional compres¬ 
sion a(t) with respect to the range of times 

t e [4000,16000], (4) 

what corresponds to hundreds of tumbling times r. We 
start late enough to avoid initial transient effects, ob¬ 
served for certain values of A. The existence of such 
transients is straightforward for fibers which were ini¬ 
tially aligned with the flow, but later remain coiled. The 
results are shown in Fig. [3] A very sharp decrease of 
Kmin and amin at Acs 0.9 defines the stiffness thresh¬ 
old which separates the coiled mode of flexible fibers from 
the straightening out mode of stiff fibers. 

Examples of the hber dynamics, discussed above for 
A = 0.7 and A = 1.14, have a common property: they 
are strictly periodic. The fiber shape at time t -I- r is 
the same as its shape at time t, rotated by the angle 
TT with respect to y-axis. Owing to this symmetry, the 
consecutive tumbling times r are equal to each other, 
and there is no net migration across the flow. Indeed, as 
shown in Figs. [Il[6),((i), a displacement Azm of the fiber 
center-of-mass position Zm at time t is compensated by 
the opposite displacement —Azm at time t+r. 

We will now analyze a different solution, called the 
migrating regular mode: superposition of periodic tum¬ 
bling with a drift across the streamlines. An example 
with A = 0.55 is shown in Figs. [T][a)-[21[a) and Movie (a). 
For such a dynamics, there is no rotational symmetry of 
shapes at time t and t + r. The consecutive tumbling 
times Ti and T 2 are not equal. The evolutions of shapes 
during ti and T 2 differ from each other, but during every 
second tumbling time are the same, with the period ti+T 2 . 
The asymmetry of the dynamics during two consecutive 
tumbling times results in a net fiber displacement along 
2 , the same after each period, as shown in Fig. [TJa). 



In Fig. m we determine the stiffness thresholds for ac¬ 
tivation of the migrating regular modes. For each value 
of A we plot all the tumbling times r within the time 
range specified in Eq. the mean migration velocity v 
(equal to the center-of-mass displacement Azm across the 
flow in the time range specified in Eq. divided by the 
time difference), and the curvature k/ = K(tf) at all the 
flipping instants tj. At Ap£ib~ 0.35 and ApppRiO.6 we 
observe period-doubling and period-halving bifurcations 
[ 2 ^ between the periodic mode (with a single value of 
the tumbling time and no migration) and the migrating 
regular mode (with two different consecutive tumbling 
times Ti and T 2 ). 

All the observed coiled solutions are periodic or migrat¬ 
ing regular, as illustrated in Eig. U) However, above the 
threshold Acs for the straightening out mode, a chaotic 
dynamics is observed. An example of such an irregular 
solution with A = 0.9 is shown in Figs. ®c), [2l[c) and 
Movie (c). A variability of the dynamics during differ¬ 
ent tumblings is visible. Subsequent tumbling times, and 
fiber curvature kj & fractional compression a/ at many 
consecutive flipping instants tf are shown in Fig. [5l Val¬ 
ues of these three quantities evolve at random between 
certain upper and lower limits, what may correspond to 
one or more unstable migrating regular solutions. The 
irregular, erratic motion of the fiber center-of-mass is il¬ 
lustrated in Fig. [6] 

Chaotic trajectories have been observed for many val¬ 
ues of the bending stiffness in the range 0.9 < A < 1.1. 
Most of them are similar as for A = 0.9, others contain 
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FIG. 3: Transition between the coiled and straightening 
out modes: the minimal fiber curvature Kmin and 
fractional compression amin versus A. 


FIG. 4: Period-doubling bifurcations between periodic 
and migrating regular modes: consecutive tumbling 
times T, migration velocity v and curvature Kf versus A. 
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FIG. 5: Irregular mode (A = 0.9): consecutive tumbling 
times T, and fiber curvature k/ & fractional 
compression aj at consecutive flipping instants tf. 


27j is equal to Ts =7r(r + l/r) =32.0. Using the hydro- 


MULTIPOLE code with the truncation at the multipole 
order L = 10, we evaluate the reference tumbling time 
Tr = 26.6 of a rigid fiber made of beads. 

Coiled flexible fibers have an effective aspect ratio 
smaller than r. Therefore, it is reasonable that their tum¬ 
bling times in Fig. |4] are below t^. However, surprisingly, 
we have found that for flexible fibers which straighten 
out, the tumbling times can exceed r^., even though their 
mean curvature is large and the mean geometrical aspect 
ratio is smaller than r. In the chaotic range, the tumbling 
times fluctuate up to r^SS, a value larger than >Trj 
and in the periodic range with A > 1.14, the tumbling 
times T>Ta>Tr. In a sense, it takes an additional time 
to change shape while tumbling. 

Concluding, the motion of flexible fibers with open 
ends and moderate aspect ratio in shear flow has been 
determined for a wide range of their bending stiffness 
0.3 < A < 50, assuming the initial alignment with the 
flow. Different modes of the dynamics have been ob¬ 
served for different values of H, as summarized in Fig. [T] 
Period doubling bifurcations and transition to chaos have 
been found, in analogy to previous studies of different 
deformable objects: flexible knotted fibers [l^ and vesi¬ 
cles dSi. 


transient parts with modulated oscillations (periodic-n 
solutions) or (occasionally) evolve towards migrating reg¬ 
ular solutions, as illustrated in Fig.[6l The evolution pat¬ 
tern is sensitive to small changes of the bending stiffness. 
In this range of A, solutions are also sensitive to the nu¬ 
merical accuracy, and therefore we do not perform any 
systematic analysis. 

For A > 1.14, solutions are periodic, with the basic 
properties discussed earlier for A = 1.14 as an example. 

It is interesting to compare the tumbling times of flex¬ 
ible fibers with two reference solutions: rigid Jeffery’s 
spheroids [131 and rigid rods made of beads {A = oo), 
both with the aspect ratios r = 10.09: the same as our 
fiber straightened out at the equilibrium position. In our 
units, the reference tumbling time Ts of the rigid spheroid 



FIG. 6: Examples of trajectories Zmit) in the chaotic 
range of the stiffness: irregular with A = 0.9 (typical) 
and migrating regular with A = 0.95 (exceptional). 



0.35 0.6 0.9 1.1-1.14 A 

FIG. 7: Dynamics of fibers with a different stiffness A. 
Modes: C - coiled, S - straightening out, p - periodic, 
m - migrating regular, ch - chaotic. Thresholds: A - 
period-doubling bifurcation, I - transition to chaos. 
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